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Abstract 

The construction of unitary operator bases in a finite-dimensional Hilbert space is re- 
viewed through a nonstandard approach combinining angular momentum theory and rep- 
resentation theory of SU(2). A single formula for the bases is obtained from a polar 
decomposition of SU (2) and analysed in terms of cyclic groups, quadratic Fourier trans- 
forms, Hadamard matrices and generalized Gauss sums. Weyl pairs, generalized Pauli 
operators and their application to the unitary group and the Pauli group naturally arise in 
this approach. 
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1 Introduction 



Angular momentum theory HI and its group-theoretical formulation in terms of the Wigner- 
Racah algebra of SU(2) [|2l[3im (see also [5J for an extension to a finite or compact group) 
are of central importance in subatomic, atomic, molecular and condensed matter physics. 
The components of any angular momentum (spin, isopin, orbital angular momentum, 
etc.) generate the Lie algebra of the group SU (2). Therefore, SU (2) and its noncompact 
extension SU(1, 1) are basic ingredients for dealing with generalized angular momenta. 
Chains of groups ending with 50(3) ~ SU(2)/Z 2 (+) or SO (3) C SO(2) are of in- 
terest in subatomic and atomic physics. In this direction, one can mention the group 
SU(3) <8> SU(2) <8> 17(1) (related to the chain 17(3) C SU(2) ® U(l) C U(l)) and its 
grand unified and/or supersymmetric extensions for describing elementary particles and 
their (strong and electroweak) interactions J6]|. Furthermore, one know the relevance in 
atomic physics of the chain U(7) C SO(7) C G 2 C 50(3) C 50(2) for the electronic 
spectroscopy of f N ions |[3). On the other side, chains ending with 577(2) C G, where 
G is a finite group (or a chain involving finite groups), proved to be of considerable in- 
terest in molecular and condensed matter spectroscopy El [SI 13. Recently, chains of type 
577(2) C G were also used in attempts to understand the flavor structure of quarks and 
leptons [[Toll . The groups SU(2) and SU(1, 1), as well as their q- or gp-deformations in 
the sense of Hopf algebras (see for instance QT1[T2]|, thus play a pivotal role in many areas 
of physical sciences. 

The representation theory of SU{2) is generally adressed in two different ways. The 
standard one amounts to diagonalise the complete set {j 2 ,j z } involving the Casimir oper- 
ator j 2 and one generator j z of SU(2). Another way is to consider a set {j 2 , v}, where v 
is an operator defined in the enveloping algebra of SU(2) and invariant under a subgroup 
of SU(2). A third way (not very well-known) consists in diagonalising a complete set 
{j 2 , v ra }, where v ra stands for a two-parameter operator which commutes with j 2 and is 
a pseudoinvariant under a cyclic group |fl3ll . 

It is the aim of this review paper to show that the third approach to the representation 
theory of SU(2) opens a window on the apparently disconnected subjects enumerated in 
the title. 

The plan of the paper is as follows. The minimal requirements for a {j 2 , v ra ] approach 
to SU(2) (i.e., a nonstandard approach to angular momentum theory) are given in Section 
2 and in two appendices. Section 3 deals with quadratic sums (in relation with quadratic 
discrete Fourier transforms, generalized Hadamard matrices, generalized quadratic Gauss 
sums and mutually unbiased bases) and Section 4 is devoted to unitary groups and Pauli 
groups. 
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The present paper is dedicated to the memory of the late Professor Yurii Fedorovich 
Smirnov who contributed to many domains of mathematical physics (e.g., Lie groups and 
Lie algebras, quantum groups, special functions) and theoretical physics (e.g., nuclear, 
atomic and molecular physics, crystal- and ligand-field theory). 

A few words about some of the notations is in order. The bar indicates complex 
conjugation. The symbol 5 a ,b stands for the Kronecker symbol of a and b. We use / and 
Id to denote the identity operator and the rf-dimensional unity matrix, respectively. The 
operator stands for the adjoint of the operator A. We note as [A, B]- and [A, B] + 
the commutator and the anticommutator of the operators A and B, respectively. We use 
the Dirac notation for a vector in an Hilbert space; furthermore, (<j)\ip) and \4>)(ip\ 
are respectively the inner product and the outer product of the vectors \ip) and \<p). The 
symbols © and stand respectively for the addition and subtraction modulo d while ® 
and l±l are used respectively for the direct product of vectors or operators and the direct 
sum of vector spaces. The matrices of type E\ jfi with the matrix elements 

( E ^)\'y := <^,a ,( W (!) 
stand for generators of the Lie group GL(d : C). For a and b coprime, we take 

+ 1 if a = k 2 mod(5) 

— 1 si a 7^ k 2 mod(6) 



(2) 



to denote the Legendre symbol of a and b (equal to 1 if a is a quadratic residu modulo b 
and —1 if a is not a quadratic residu modulo b). In addition, the integer inverse (a\b) of a 
with respect to b is given by 

a(a\b) = 1 mod(6). (3) 

Finally, the g-deformed number [n] and the g-deformed factorial [n] \, with n E N, are 
defined by 

1 - a n 

and 

[n\ q \ := [l\ q [2\ q . . . [n\ q [0] ff ! := 1 (5) 
where q is taken is this paper as a primitive root of unity. 
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2 A nonstandard approach to su(2) 



In some previous works [fT3l . we developed a nonstandard approach to the Lie algebra 
su(2) and studied the corresponding Wigner-Racah algebra of the group SU(2). This 
nonstandard approach is based on a polar decomposition of su(2), based in turn on a 
troncated oscillator algebra (see Appendices A and B). It yields nonstandard bases for the 
irreducible representations of SU (2) and new Clebsch-Gordan coefficients for the angular 
momentum theory. Basically, the approach amounts to replace the set {j 2 , j z }, familiar in 
quantum mechanics, by a set {j 2 , v ra ] (j 2 and j z are the Casimir operator and the Cartan 
generator of su(2), respectively). 

The operator v ra acts on the (2j + 1) -dimensional subspace £(2j + 1), associated with 
the angular momentum j, of the representation space of SU{2). We define it here by 

3-1 

Vra := e^ r \j, + £ q ( °- m)a \3,m+ l)(j,m\ (6) 

m=-j 

where 

g:=eXP (27Tl) 2jeN aGZ 2,+i (7) 

and, for fixed j, the vectors \j, m) (with m = j, j — 1, . . . , — j) satisfy the eigenvalue 
equations 

j 2 \j } m) = j(j + l)\j,m) j z \j,rn) = m\j,m) (8) 

familiar in angular momentum theory. The vectors \j, m) span the Hilbert space £(2j + 
1) ~ C 2j+1 and are taken in an orthonormalized form with 

(j,m\j,m') = 5 m>m i. (9) 

Obviously, the operator v ra is unitary and commutes with j 2 . The spectrum of the set 
{j 2 ! v ra } is described by 

Result 1. For fixed j, r and a, the 2j + 1 vectors 

1 j 

\ja;ra) := , ^ q^ +m ^- m+l ^ 2 -^ mr+ ^ +m ^ a \j,m) (10) 

with a = 0, 1, ... , 2j, are common eigenvectors ofv ra and j 2 . The eigenvalues ofv ra and 
j 2 are given by 

v ra \ja;ra) = q j{a+r) ~ a \ja;ra) j 2 \ja;ra) = j(j + l)\ja;ra) a = 0, 1, . . . , 2j. (11) 
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The spectrum ofv ra is nondegenerate. 

The set {\ja; ra) : a = 0, 1, . . . , 2j} constitutes another orthonormal basis, besides 
the basis {\j, m) : m — j, j — 1, . . . , — j}, of S(2j + 1) in view of 

(J a; ra\jft; ra) = S a)j3 . (12) 

Note that the value of (ja; ra\j{3; sb) is much more involved for r ^ s and a ^ b and 
needs the calculation of Gauss sums as we shall see below. 

The Wigner-Racah algebra of SU (2) can be developed in the {j 2 , v ra ] scheme. This 
leads to Clebsch-Gordan coefficients and (3 - ja) ra symbols with properties very differ- 
ent from the ones of the usual SU(2) C U(l) Clebsch-Gordan coefficients and 3 - jm 
symbols corresponding to the {j 2 ,j z } scheme Ifl3"l . 

The nonstandard approach to angular momentum theory briefly summarized above is 
especially useful in quantum chemistry for problems involving cyclic symmetry. This is 
the case for a ring- shape molecule with 2j + 1 atoms at the vertices of a regular poly- 
gon with 2j + 1 sides or for a one-dimensional chain of 2j + 1 spins (|-spin each) |[T4l . 
In this connection, we observe that the vectors of type \ja; ra) are specific symmetry- 
adapted vectors [|T5l[T6l . Symmetry-adapted vectors are widely used in quantum chem- 
istry, molecular physics and condensed matter physics as for instance in rotational spec- 
troscopy of molecules ifTTl and ligand-field theory lfT8l . However, the vectors \ja;ra) 
differ from the symmetry-adapted vectors considered in Refs. lfT9ll20ll2~Tll22ll in the sense 
that v ra is not an invariant under some finite subgroup (of crystallographic interest) of the 
orthogonal group 0(3). Indeed, v ra is a pseudoinvariant [|23l under the Wigner operator 
PrUp) associated with the rotation R(ip), around the quantization axis Oz, with the angle 

ip:= P7 ^- p = 0,l,...,2j (13) 
2j + 1 

since 

P%)^ra4 W = e~^V ra . (14) 

More precisely, we have 

Result 2. The operator v ra transforms according to an irreducible representation of the 
cyclic subgroup C^j+i ~ Z 2j+ i(+) of the special orthogonal group SO (3). In terms of 
vectors, one has 

■Pjafoob'a; ra) = q w \jp; ra) f3 := a Qp (15) 

so that the set { \ja; ra) : a — 0, 1, . . . , 2j} is stable under Pr^). The latter set spans the 
regular representation ofC 2 j+i- 
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3 Variations on quadratic sums 
3.1 Quadratic discrete Fourier transform 

We leave the domain of angular momentum theory and adopt the following notations 

d := 2j + 1 k:=j-m \k) := \j,m). (16) 

These notations are particularly adapted to quantum information and quantum computa- 
tion. In these new notations, we have 

d-l 

V ra = e^-^ld - 1)(0| + Y, ( l ka \ k - l )( k V ( 17 ) 

k=l 

From now on, we assume that d > 2 and r = (the case d — 1 and r / 0, although of 
interest in the theory of angular momentum, is not essential for what follows). In addition, 
we put 

\aa) := \ja; 0a) (18) 
with a and a in the ring Z d := Z/dZ. Then, Eq. (fTTT) gives 

v 0a \aa) = q {d - 1)a/2 - a \aa) (19) 

with 

| aa ) = _L ^^(fc+iJCd-fc-Da/a-Cfc+Da^^ (20) 
Equation (|2Q|) can be rewritten as 

«> = £( F «)*«I*> (2D 



c/n ) 

fc=0 



where 

/p\ ._ 1 g (fc+i)(d-fc-i)o/2-(fc+i)a ^2) 

is the /ca-th matrix element of a d x d matrix F a (the matrix F a can be seen as the matrix 
associated with the character table of the cyclic group Cd pre- and post-multiplied by 
diagonal matrices). 
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Equations (|2TI)-(|221) define a quadratic quantum Fourier transform. The matrix F a is 
unitary so that (|2TT) can be inverted to give 



d-l 



\k} = Y,( F ^kJ^) (23) 



a=0 

or 

\k) = 4= Vg- (fc+1)(d " fc - lW2+(fe+1)a |a«). (24) 
In the special case a = 0, we have 

|0a) = g" Q 4 52 ^ |jfc) = 4 5" e^ fc+1 > a |0a>. (25) 

v fe=0 a=0 

Consequently, the quadratic quantum Fourier transform reduces to the ordinary quantum 
Fourier transform (up to a phase factor). The corresponding matrix F satisfies 

F% = qh (26) 

to be compared to the well-known relation F 4 = I d for the standard quantum Fourier 
transform [|24l . 

At this stage, we forsee that d + 1 (orthonormal) bases of the space £{d) play an 
important role in the present paper: (i) the basis 

B d := {\j, m):m = j,j-l,...,-j}&B d :={\k):k = 0,l,...,d-l} (27) 

associated with the {j 2 ,j z } scheme, known as the spherical or canonical basis in the 
theory of angular momentum, and as the computational basis in quantum information and 
quantum computation and (ii) the d bases 

B a := {|aot) : a = 0,1,.. .,d- 1} a = 0, 1, . . . , d - 1 (28) 

(noted B 0a in Ref. 11251 ) associated with the {j 2 , v 0a } scheme. 

To close this subsection, let us show how the preceding developments can be used for 
defining a quadratic discrete Fourier transform. We start from the formal transformation 

x ■= { x (k) e C : k = 0, 1, . . . , d - 1} -»• y := {y(a) e C : a = 0, 1, . . . , d - 1} (29) 

defined via 

d-l 

Vd 



y ( a ) := ^ S g (fc+1)(d - fc - 1)a/2 - (fe+1)a x(A;) (30) 



fc=0 
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where a can take any of the values 0, 1, . . . , d — 1. Alternatively, for fixed a we have 

y(a) = (F a ) ka x(k) a = 0,1,..., d-l. (31) 

fc=0 

The inverse transformation y — > x is described by 

d-l 

<k) = YjFahAa) k = 0,1,..., d-l. (32) 

o=0 

The bijective transformation x <-> y can be thought of as a quadratic discrete Fourier 
transform. The case a = corresponds to the ordinary discrete Fourier transform (up to a 
phase factor). These matters lead to the following result which generalizes the Parseval- 
Plancherel theorem for the ordinary discrete Fourier transform. 

Result 3. The quadratic discrete Fourier transforms x *-> y and x' <-> y' associated withe 
same matrix matrix F a , a G Z d , satisfy the conservation rule 

d-l d-l 

Y^)v\oc) = Y^x'ik) (33) 

o=0 fc=0 

where the common value is independent of a. 

3.2 Generalized Hadamard matrices 

The modulus of each matrix element of F a (with a G Z^) is equal to 1/ \fd. Therefore, 
the unitary matrix F a turns out to be a generalized Hadamard matrix. We adopt here the 
following definition. A d x d generalized Hadamard matrix is a unitary matrix whose 
each entry has a modulus equal to 1/y/d G6l . Note that the latter normalization, used 
in quantum information If27ll28l , differs from the usual one according to which a d x d 
generalized Hadamard matrix if is a complex matrix such that H^H = did and for which 
the modulus of each element is 1 ||29l . In this respect, the generalized Hadamard matrix 
H a considered in lfl4l corresponds to ydF a up to permutations. 

Example 1. By way of illustration, from (|22|) we get the familiar Hadamard matrices 



f i = — [ ; j I (34) 



for d = 2 and 





to 2 ' 




1 J F 2 = — J ,c 2 1 a: J (35) 
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(with uj : = e* 27r//3 ) for d = 3. Another example is 

-2 



/ 1 T T 

2 



F ° = 7E 



1 -r 

.2 



r 2 \ 

1 r* -r 1 r 2 -r 
1-1 1-11 -1 
1 -r r 2 1 -r r 2 
2 -r -1 r 2 r 
111 1 / 



(36) 



1 — T' 
V 1 1 

(with r := e~ l7r//3 ) which readily follows from (1221) for d = 6 and a = 0. 



We sum up and complete this section with the following result (see also lTl4l[30ll '). 
Result 4. The matrix 

F a = -L E g( fc+1 )( d - fc - 1 ) fl / 2 -( fc+1 ) Q E fc , a (37) 

V « fc=0 «=0 

associated with the quadratic quantum Fourier transform ri271) is a d x d generalized 
Hadamard matrix. It reduces the endomorphism associated with the operator t> 0a -' 



d-l 



FlV 0a F a = q^ 2 ^~ a Ea, a = q {d - 1)a 
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a=0 



(\ 

q- 1 

\0 



\ 


-(d-l) 



(38) 



where the matrix 



d-l 



feei.fc 



fc=0 



/0 q a 

q 2a 



\1 



Q 



\ 


(d-l)a 

o y 



(39) 



represents the linear operator v 0a on the basis B d . 



3.3 Generalized quadratic Gauss sums 

The Hadamard matrices F a and F b (a, b G Z d ) are connected to the inner product (aa\b(3). 
In fact, we have 



(aa\bf3) = (Fift) 



A direct calculation yields 



( aa |6/3) = 



fc(d-ifc)(b-o)/2-fe(/3-a) 



(40) 



(41) 



fc=0 



or 

d-l 

(aa\b/3) = - ^ e M(»-6)fe 2 +[d(6-a)+2(a-/3)]fe}/d_ ^ 
^ fc=0 

Hence, each matrix element of F^F b can be put in the form of a generalized quadratic 
Gauss sum S(u, v, w) defined by QTTl 

\w\-l 

S(u,v,w) := e m{uk2+vk)/w (43) 

k=0 

where u, v and w are integers such that u and w are mutually prime, uw ^ and uw + v 
is even. In detail, we obtain 

(aa\b{3) = {FlF b ) ap = - d S{u 1 v ) w) (44) 

with the parameters 

u = a — b v = —(a — b)d + 2(a — (3) w = d (45) 

which ensure that uw + v is necessarily even. 

In the particular case d = 2 (of special interest for qubits), we directly get 

(aa\b(3) = -[! + e ^( b -+^-^] (46) 

which reduces to 

(aa\a(3) = 5 a ^ b = a, (47) 

and 

(aa\b/3) = ^(l±i) b ^ a, (48) 

where the + sign corresponds to 6 — a + 2(a — (3) = 1,-3 and the — sign to b — a + 
2{a-(3) = -1,3. 

In the general case d arbitrary (of interest for qudits), the calculation of S(u, v, w) can 
be achieved by using the methods described in PTfl (see also (321 [331 ES [351). The cases 
of interest for what follows are (u even, v even, w odd), (u odd, v odd, w odd) and (u odd, 
v even, w even). This leads to 

Result 5. For a ^ b, d arbitrary and u + v + w odd, the inner product {aa\bf3) and the 
a(3-th element of the matrix F^F b follow from 
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case u = a — b even, v = d(b — a) + 2(a — 0) even, w = d odd: 




case u = a — b odd, v = d{b — a) + 2 (a — (3) odd, w = d odd: 




case u = a — b odd, v = d(b — a) + 2 (a — (3) even, w = d even: 

(aa\b(3) = (FiF b ) ap = ^ (^exp + l(n\ W )V]) (51) 

50 matrix F^Ff, is a Hadamard matrix for each case under consideration. 

Finally, for a = b and d arbitrary we recover the orthonormality property (see (fT2~l) ) 

(aa\a/3) = 5 a ,p (52) 
from a direct calculation of the right-hand side of (l42j) . 

3.4 Mutually unbiased bases 

Speaking generally, two <i-dimensional bases B a = {\aa) : a G Z d } and B b = {\b(3) : 
(3 G Z d } are said to be mutually unbiased if and only if 

\(aa\bp)\ = 6 atb 6 a ,p + (1 - (53) 

for any a and /3 in the ring 7h d . It is well-known that the number of mutually unbiased 
bases (MUBs) in the Hilbert space C d cannot be greater than d + 1 Ii36l l37l l38l [39ll . In 
fact, the maximum number d + 1 is attained when d is the power of a prime number 
j|38l l39ll . Despite a considerable amount of works, the maximum number of MUBs is 
unknown when d is not a power of a prime. In this respect, several numerical studies 
strongly suggest that there are only three MUBs for d = 6 (see for example E71 l28l l40l 
SHI421). MUBs are closely connected with the concept of complementarity in quantum 
mechanics. There are of paramount importance in classical information theory (Kerdock 
codes and network communication protocols) [|39l |43l , in quantum information theory 
(quantum cryptography and quantum state tomography) [|44|| and in the solution of the 
Mean King problem [@5l|46l|47l|48l|49l[50l|. Recently, it was pointed out and confirmed 
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that MUBs are also of central importance in the formalism of Feynman path integrals 
||5Tl[52l . Finally, it should be emphasized that the concept of MUBs also exists in infinite 
dimension 11531 . There are numerous ways of constructing sets of MUBs. Most of them 
are based on discrete Fourier analysis over Galois fields and Galois rings, discrete Wigner 
functions, generalized Pauli matrices, mutually orthogonal Latin squares, finite geometry 
methods and Lie-like approaches (see Refs. 01|21llSl|2SllSaE3E2lESlE5lffiaBIl 

H3 SH SE1 SSI El ISDl H31 US Hzl HH1 HH IgOl [gH IS1 for an 

nonexhaustive list of references). 
3.4.1 Case d prime 

We have the following important result. (See also [64] for a recent alternative group- 
theoretical approach to the case d prime.) 

Result 6. In the case where d = p is a prime number ( even or odd), one has 



for a, b,a, (3 G Z p . Therefore, the p + 1 bases B , Bi,...,B p constitute a complete set of 
MUBs in CP. 

The proof easily follows from the calculation of the modulus of S(a — b,pb — pa + 
2a - 2(3, p) from (|46b . ( 149b , (l50b and (|5_D. As a consequence, the bases B a , with a = 
0, 1, . . . , p — 1, are p MUBs in the sense that they satisfy (l53b for any a, b, a and (3 in the 
Galois field F p . Obviously, each of the bases B a (with a — 0, 1, ... ,p — 1) is mutually 
unbiased with the computational basis B p . This completes the proof. Note that Result 6 
can be proved as well from the developments in [fT4l . 

As two typical examples, let us examine the cases d = 2 and 3. 

Example 2: case d = 2. In this case, relevant for a spin j = 1/2 or for a qubit, we have 
q = — 1 and a, a E Z 2 . The matrices of the operators v 0a are 



\(aa\b(3)\= (F*F b ) 



1 



a b 



(54) 




(55) 



By using the notation 




1 1 



(56) 
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familiar in quantum chemistry (a is a spinorbital for spin up and (3 for spin down), the 
d + 1 = 3 MUBs are 

B : \00) = ^=(a + P) |01) = ~ (a - (3) (57) 

B 1 : \10) = (a - i/3) \U) = -i^= (a + i/3) (58) 
v 2 v 2 

5 2 : |0) = a |l)=/3. (59) 



Example 3: case d = 3. This case corresponds to a spin j = 1 or to a qutrit. Here, we 
have g = exp(i27r/3) and a, a G Z 3 . The matrices of the operators t> 0a are 



oo 




Vr 



01 



The d + 1 = 4 MUBs read 




Vr 



02 




(60) 



Bn 



B 1 



Bo 



00) 
01) 
02) 
10) 

11) 
12) 

20) 

21) 

22) 
0) = 



^(|0> + |1> + |2» 
i=(g 2 |0) + g|l) + |2)) 



(g|0) + g 2 |l) + |2)) 



_(g|0>+ (Z |l> + |2» 

-L(|0) + g 2 |l) + |2)) 



(g 2 |0) + |l) + |2)) 



- 7 =(g 2 |0) + g 2 |l) + |2)) 



1 

i 

"7s 
IM> 



(?|0> + |1> + |2» 

(|0> + «z|l> + |2» 

|1) = |1,0) |2> = |1,-1>. 



(61) 



(62) 



(63) 
(64) 



It should be observed that B (respectively, B\ and B 2 ) can be associated with the vector 
(respectively, projective) irreducible representations of the group C 3 . 
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3.4.2 Case d power of a prime 

Different constructions of MUBs in the case where d is a power of a prime were achieved 
by numerous authors from algebraical and geometrical techniques (see for instance ll38l 
E9lllllS3ll3ll2llSHl|55l|gDll5Il|g2llSaand references therein). We want to show here, 
through an example for d = 4, how our angular momentum approach can be useful for 
addressing this case. 

Example 4: case d — 4. This case corresponds to a spin j = 3/2. Here, we have q = i 
and a, a E Z 4 . Equations (|20l) and (|28l ) can be applied to this case too. However, the 
resulting bases B , B\, B 2 , B 3 and £> 4 do not constitute a complete system of MUBs 
(d = 4 is not a prime number). Nevertheless, it is possible to find d + 1 = 5 MUBs 
because d = 2 2 is the power of a prime number. This can be achieved by replacing the 
space £(4) spanned by {|3/2, m) : m = 3/2, 1/2, —1/2, —3/2} by the tensor product 
space 8(2) <g> £(2) spanned by the basis 

{a®a,a®l3,f3®a,f3®f3}. (65) 

The space 8(2) ® 8(2) is associated with the coupling of two spin angular momenta 
ji = 1/2 and j% = 1/2 or two qubits (in the vector u <g> v, u and v correspond to j\ and 
j 2 , respectively). 

In addition to the basis (|65l) , it is possible to find other bases of 8(2) Cg> 8(2) which are 
mutually unbiased. The d = 4 MUBs besides the canonical or computational basis (|65l) 
can be constructed from the eigenvectors 

\aba(3) := \aa) ® \bf3) (66) 

of the operators 

Wab ■= V 0a ® Wof, (67) 

(the vectors |aa) and |&/3) refer to the two spaces 8(2)). As a result, we have the d+ 1 = 5 
following MUBs where A = (1 — z)/2 and /i = i\. 
The canonical basis: 

a® a a®f3 j3 ® a f3®f3. (68) 

The woo basis: 

|0000) = \(a®a + a®p + P®a + fi®P) (69) 
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(82) 


A|1001> + /x|1010> = 
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g)/3 + /3< 


8) a 




5/3) 


(83) 


^1001) + A|1010) = 




8 a 


— ia ( 


8/5 


-/3 


8 a 




5/3). 


(84) 



It is to be noted that the vectors of the w o and wu bases are not intricated (i.e., each 
vector is the direct product of two vectors) while the vectors of the w i and w w bases are 
intricated (i.e., each vector is not the direct product of two vectors). To be more precise, 
the degree of intrication of the state vectors for the bases w 00 , w n , w 01 and w 10 can be 
determined in the following way. In arbitrary dimension d, let 

d-l d-l 

|$) = ^$> fe /|A:)®IO (85) 

k=0 1=0 
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be a double qudit state vector. Then, it can be shown that the determinant of the d x d 
matrix A = (a w ) satisfies 

< IdetAI < --L (86) 

" Vd 1 



as proved in the Albouy thesis 11651 16611 . The case det A = corresponds to the absence 
of intrication while the case 

I det A| = — = (87) 
Vd* 

corresponds to a maximal intrication. As an illustration, we obtain that all the state vec- 
tors for w o and wn are not intricated and that all the state vectors for u> i and wio are 
maximally intricated. 

3.4.3 Case d arbitrary 

In the special case where u — 1, the generalized Gauss sum 5(1, —d + 2a — 2/3, d) can 
be easily calculated for d arbitrary by means of the reciprocity theorem [[311 

S(u, V,w) = \ n^ e Msgz(uw)-v*/(uw)]/4 S f Wj u) _ (8g) 

V u 

This leads to the following particular result. 

Result 7. For d arbitrary and b = a 1, one has 

(aa|a_!/3) = ^ e ™[i-(d-2 a+ 2pf /j\/a ^ | ( aa | a _ l/3 ) | = J_ a _i = o e 1. (89) 
yd yd 

Therefore, the three bases B ae i, B a and are mutually unbiased in C d . 

This result is in agreement with a well-known result proved in many papers from 
quite distinct ways (see for instance PTi ). We thus recover, from an approach based on 
generalized Gauss sums, that for d arbitrary the minimum number of MUBs is 3. 

Another special case, viz., u = 2 (=> d > 3), is worth of value. The application of the 
reciprocity theorem gives here 

Result 8. For d > 3 and b = a 2, one has 

(aa\a_ 2 (3) = ±= -1= ^-^-^/^[i + e <7r(-d+a«-2/j)/2i 
VdV2 



\{aa\a-2p)\ = \l — cos — (d — 2a + 2(3) 



a„ 2 = a 2.(90) 
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Therefore, the bases B aQ 2 and B a cannot be mutually unbiased in C d for d even with 
d > 4. In marked contrast, the bases B ae2 and B a are unbiased for d odd with d > 3 ( d 
prime or not prime). 

Going back to the Hadamard matrices, let us remark that, for d arbitrary, if B a and 
B b are two MUBs associated with the Hadamard matrices F a and F b (respectively), then 
F^Fb is a Hadamard matrix too. However, for d arbitrary, if F a and Ft, are two Hadamard 
matrices associated with the bases B a and B b (respectively), the product F\F b is not in 
general a Hadamard matrix. 

4 Unitary group and generalized Pauli group 
4.1 Weyl pairs 

We continue with the general case where d is arbitrary. The operator v 0a can be expressed 
as 

d-l j 

v 0a = J2<l ka \kel)(k\^v 0a = q U ' m)a \j,m®l)(j,m\ (91) 

k=0 m=—j 

so that 

v 0a \k) = q ka \kQl)^ v 0a \j, m) = q^ a \j, m ® 1) (92) 

where q = exp(2ni/d). The operators x (the flip or shift operator) and z (the clock 
operator), used in quantum information and quantum computation (see for instance ||67l 
68]), can be derived from the generic operator v Qa as follows 

x:=v 00 z := (u o) t uoi- (93) 

Therefore, we get 

x = J2\kQl)(k\ = |d-l)(0| + |0)(l| + ... + |d-2)(d-l| (94) 

k=0 

and 

d-l 

z = Y^q k \k)(k\ = |0)(0|+g|l)(l| + ... + g d - 1 |c/-l)(^-l|- (95) 

k=0 

The action of x and z on the basis B d of S(d) is given by the ladder relation 

x\k) = \k e 1) <£> x\j,m) = (1 - 5 mJ ) \j,m+ 1) + 5 m>j \j, -j) (96) 
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and the phase relation 

z\k) = q k \k) z\j, m) = q j - m \j, m). (97) 
Alternatively, the action of x and z on any basis B a (a = 0, 1, . . . , d — 1) of £{d) reads 
x\aa) = q [d ' 1)a/2 ~ a \aa a ) a a = a © a x\0a) = q~ a \0a) (98) 

and 

z\aa) = g _1 |aa_i) a_i = a 1. (99) 



Equations (1961 ) and (1971 ), on one side, and Eqs. (1981) and (1991) , on the other side, show that 
the flip or clock character for x and z is basis-dependent. The relationship between x and 
z can be understood via the following 



Result 9. The unitary operators x and z are cyclic and q-commute: 



They are connected by 



x d = z d = I xz — qzx = 0. 



x = pzf <^ z = fxf* 



(100) 



(101) 



where the Fourier operator 



1 d-l d-1 



(102) 



is unitary and satisfies 



k=0 k'=0 



/* = !• 



r/ze operators x and z are isospectral operators with the common spectrum {1 , q, . 
A direct proof of Result 9 can be obtained by switching to the matrices 



(103) 

,q d - 1 }- 





(o 


i 


. 


• °\ 


d-l 








1 . 


. 












k=Q 
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. 1 









. 


• 0/ 



(104) 
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d-1 



q k Ek,k 



/I 







Q 















(105) 



fc=0 



\0 



q 



V 



of the operators x and z, in the basis B d (cf. (1971 ) and (1981)). Let F be the matrix of the 
linear operator / in the basis B d . The reduction by means of F of the endomorphism 
associated with the matrix X yields the matrix Z. In other words, the diagonalization of 
X can be achieved with the help of the matrix F via Z = FXF^. Note that the matrix F 
is connected to F by 



S :-- 



0=0 



(106) 



where S acts as a pseudopermutation. 

In view of (HOOK the pair (x, z) is called a Weyl pair. Weyl pairs were originally 
introduced in finite quantum mechanics [|69l and used for the construction of unitary bases 
in finite-dimensional Hilbert spaces irTOl . It should be noted that matrices of type X and 
Z were introduced long time ago by Sylvester [fTTI in order to solve the matrix equation 
PX = XQ; in addition, such matrices were used by Morris 11721 to define generalized 
Clifford algebras in connection with quaternion algebras and division rings. Besides the 
Weyl pair (x,z), other pairs can be formed with the operators v 0a and z. Indeed, any 
operator v 0a (a G Z d ) can be generated from x and z since 



v 0a = xz . 

Thus, Eq. (1 1001 ) can be generalized as 

e -Md-i)a( VQa y = z d = I VQaZ _ qzVQa = 

Therefore, the pair (v 0a , z) is a Weyl pair for (d — l)a even. 



(107) 



(108) 



4.2 Generalized Pauli matrices 

For d — 2 the g-commutation relation of x and z reduces to an anticommutation relation. 
In fact, Eq. (1 1001 ) with d = 2 can be particularized to the relations 

x 2 = z 2 = I xz + zx = (109) 

which are reminiscent of relations satisfied by the Pauli matrices. Hence, we understand 
that the matrices X and Z for d arbitrary can be used as an integrity basis for producing 
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generalized Pauli matrices 03 11HI IIS 113 EZl I1H1 Illl ESI Ell ISl Ell 1771 ESI 123 
[801 IBTl [82ll . Let us develop this point. 
For d arbitrary, we define the operators 

u ab = x a z h a,beZ d . (110) 

The operators u ab shall be referred as generalized Pauli operators and their matrices as 
generalized Pauli matrices. They satisfy the ladder-phase relation 

Uab\k) = q kb \kea) ^u ab \j,m) = q {j - m)b \j,m © a) (111) 

from which we can derive the following result. 

Result 10. The d 2 operators u ab , with a, b G Z^, are unitary and obey the multiplication 
rule 

UabUa'b' = U a "b" 0>" '■= O ® O b" \= b © V ' . (H2) 

Therefore, the commutator and the anticommutator ofu ab and u a /y are given by 

[u ab , u a , v \± = (q- ba ' ± q- aV ) u a „y, a" := a © a' b" := b © b' . (113) 
Furthermore, they are orthogonal with respect to the Hilbert-Schmidt inner product 

Tr £ ( d ) [{u ab ) ] u a > b ] =d8 a , a 'h,v (114) 
where the trace is taken on the d-dimensional space £(d). 
As a corollary of Result 10, we have 

[u ab ,u a / b /]- = <=>■ ab' bo! = (115) 

and 

[uab, u a 'b>] + = <^> ab' © bo! = ^d. (116) 

This yields two consequences. First, Eq. (II 161) shows that all anticommutators [u ab , u a ' b '} + 
are different from if d is an odd integer. Second, from Eq. (II 151 ) we have the important 
result that, for d arbitrary, each of the three disjoint sets 



e . 


:= {u 0a = z a : a 


= 1,2,.. 


,d-l} 


(117) 


e.. 


■ \u aa X z . 


a = 1,2, 


...,d-l} 


(118) 


e.o 


:= {u a0 = x a : a 


= 1,2,.. 


,d-l} 


(119) 
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consist of d — 1 mutually commuting operators. The three sets e „ e„ and e. are asso- 
ciated with three MUBs. This is in agreement with the fact that the bases B , B x and B d 
are three MUBs for d arbitrary (v 00 = x G e #0 , v 01 = xz G e„ and z G e . are associated 
with £? , i?x and B d , respectively). 

By way of illustration, let us give the matrices in the basis B d of the operators u ab for 
d = 2, 3 and 4. 

Example 5: case d = 2. For d = 2 <^ j = 1/2 (=>• g = —1), the matrices in the two sets 

£ := {I 2 = X°Z ,X = X 1 Z° = V 00 } (120) 
: = {Z = X°Z\Y = X 1 Z 1 = V 01 } (121) 

corresponding to the four operators u a b are 

;) i) «-(; -0 "-(: o 1 ) 

In terms of the usual (Hermitian and unitary) Pauli matrices a x , a y and a z , we have 

X = a x Y = —ic y Z = a z . (123) 

The matrices X, Y and Z are thus identical to the Pauli matrices up to a phase factor for 
Y. This phase factor is the price one has to pay in order to get a systematic generalization 
of Pauli matrices in arbitrary dimension. 

Example 6: case d = 3. For d = 3 j = 1 (=>- g = exp(i27r/3)), the matrices in the 
three sets 

(124) 
(125) 
(126) 

corresponding to the nine operators u ab are 

X = | (i I I A"- = j I \ (127) 



i? 


:= {X ^ ^ 1 ^ 


= V 00 ,X 2 Z } 




:= {X*Z X ,X X Z X 


= V 01 ,X 2 Z 1 } 




:= {X°Z 2 : X 1 Z 2 


= V 2,X 2 Z 2 } 






1 











Q 











Q 2 



XZ 
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g 2 

1 



X 2 Z 



g 2 ' 

1 
q 



(128) 
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xz 1 




2 ry2 



x A z 




■ (129) 



These generalized Pauli matrices differ from the Gell-Mann matrices and Okubo matrices 
used for SU (3) in particle physics with three flavors of quarks (831 IHH [EH. They consti- 
tute a natural extension based on Weyl pairs of the Pauli matrices in dimension d = 3. 

Example 7: case d — 4. For d = 4 j = 3/2 (=>■ q — i), the matrices in the four sets 



£7 
Si 

S3 



{X^X 1 ^ = V 00 ,X 2 Z°,X 3 Z } 
{X°Z\X 1 Z 1 = V 01 ,X 2 Z 1 ,X 3 Z 1 } 

{x°z 2 ,x l z 2 = K)2,x 2 z 2 ,x 3 z 2 } 

{X ^ 3 ^ 1 ^ 3 ee 1/ 03 , X 2 Z 3 , X 3 Z 3 } 



(130) 
(131) 
(132) 
(133) 



corresponding to the 16 operators u ab are 
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(134) 



(135) 



(136) 



X 2 Z 



/0 -1 \ 

-i 

10 

\0 i 1 / 



X 6 Z 



/0 -A 

10 

i 

\0 -1 / 



(137) 



/l \ 
0-100 

10 
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(138) 
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2 ryl 
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(139) 



(140) 
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(141) 



These generalized Pauli matrices are linear combinations of the generators of the chain 
S77(4) d 5C/(3) D 5C/(2) in particle physics with four flavors of quarks ||86l[87l[88l. 



For d arbitrary, the generalized Pauli matrices arising from (111 Ob are different from 
the generalized Gell-Mann A matrices introduced in [89]. The generalized A matrices are 
Hermitian and adapted to the chain of groups SU(d) D SU(d — 1) D . . . D SU(2) while 
the matrices X a Z b are unitary and closely connected to cyclic symmetry. Indeed, for d 
arbitrary, each of the d sets 



E h 



{X a Z b 



0,1, 



d-1} 



0,1, 



d-1 



(142) 



is associated with an irreducible representation of the cyclic group Ca- More precisely, 
the one-dimensional irreducible representation of Cd associated with Eb is obtained by 
listing the nonzero matrix elements of any matrix of E b , column by column from left to 
right. In this way, we obtain the d irreducible representations of Cd- This relationship 
between rf-dimensional Pauli matrices and irreducible representations of Cd are clearly 
emphasized by the examples given above for d = 2, 3 and 4. 



4.3 Pauli basis for the unitary group 



Two consequences follow from (II 141) . (i) The Hilbert-Schmidt relation (II 141) in the 
Hilbert space C rf2 shows that the d 2 operators u a b are pairwise orthogonal operators. Thus, 
they can serve as a basis for developing any operator acting on S(d). (ii) The commuta- 
tor in (II 131) defines the Lie bracket of a d 2 -dimensional Lie algebra generated by the set 
{u ab : a, b — 0, 1, . . . , d— 1}. This algebra can be identified to the Lie algebra u(d) of the 
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unitary group U(d). The subset {u a b : a,b — 0,1, ... ,d — 1} \ {u 00 } then spans the Lie 
algebra su(d) of the special unitary group SU (d). In other words, the Weyl pair (X, Z), 
consisting of the generalized Pauli matrices X and Z in dimension d, form an integrity 
basis for u(d). More specifically, the two following results hold. 

Result 11. The set {X a Z b : a, b = 0, 1, . . . , d — 1} forms a basis for the Lie algebra u{d) 
of the unitary group U (d). The Lie brackets ofu(d) in such a basis (denoted as the Pauli 
basis) are 

d-l d-l 

[x «z b ,X e Zf}_ = ^^(a&,e/;u)^' (143) 

with the structure constants 

(ab, ef; ij) = 5(i, a © e)S(j, b © /) (q- be - q'^) (144) 

where a, b, e, f, i, j G Z d . The structure constants (ab, ef; ij) with i = a©e and j = b®f 
are cyclotomic polynomials associated with d. They vanish for af be = 0. 

Result 12. For d — p, with p a prime integer, the Lie algebra su (p) of the special 
unitary group SU (p) can be decomposed into a direct sum ofp+1 abelian subalgebras 
of dimension p — 1, i.e. 



su(p) ~ v Q l±l vi l±l . . . I±l v p (145) 

where each of the p+1 subalgebras Vq, v±, . . . , v p is a Cartan subalgebra generated by a 
set ofp — 1 commuting matrices. The various sets are 

V := {X°Z 1 ,X Z 2 ,X Z 3 ,...,X Z p - 2 ,X°Z p - 1 } (146) 

Vi := {X 1 Z°,X 2 Z ,X 3 Z ,...,X p - 2 Z°,X p - 1 Z } (HI) 

V 2 := {X 1 Z 1 ,X 2 Z 2 ,X 3 Z 3 ,...,X p - 2 Z p - 2 ,X p - 1 Z p - 1 } (148) 

V 3 := {X 1 Z 2 ,X 2 Z\X 3 Z 6 ,...,X p - 2 Z p - 4 ,X p ' 1 Z p - 2 } (149) 

; (150) 

V p _i := {X 1 Z p - 2 ,X 2 Z p -\X 3 Z p - 6 ,...,X p - 2 Z\X p - 1 Z 2 } (151) 

V p := {X 1 Z p - 1 ,X 2 Z p - 2 ,X 3 Z p - 3 ,...,X p - 2 Z 2 ,X p - 1 Z 1 } (152) 



for Vq, Vi, . . . , v p , respectively. 
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Example 8: p = 7 <^> j = 3. Equations (I146I) - (I152I) give 



V = {(01), (02), (03), (04), (05), (06)} (153) 

Vi = {(10), (20), (30), (40), (50), (60)} (154) 

V 2 = {(11), (22), (33), (44), (55), (66)} (155) 

V 3 = {(12), (24), (36), (41), (53), (65)} (156) 

V 4 = {(13), (26), (32), (45), (51), (64)} (157) 

V 5 = {(14), (21), (35), (42), (56), (63)} (158) 

V 6 = {(15), (23), (31), (46), (54), (62)} (159) 

V 7 = {(16), (25), (34), (43), (52), (61)} (160) 

where (ab) is used as an abbreviation of X a Z b . 



Result 12 can be extended to the case where d = p e with p sl prime integer and e a 
positive integer: there exists a decomposition of su(p e ) into p e + 1 abelian subalgebras of 
dimension p e — 1. In order to make this point clear, we start with a counterexample. 

Counterexample: d = 4 4» j — 3/2 (=> a, 6 = 0, 1, 2, 3). In this case, Result 1 1 is valid 
but Result 12 does not apply. Indeed, the 16 unitary operators u a b corresponding to 

ab = 01, 02, 03, 10, 20, 30, 11, 22, 33, 12, 13, 21, 23, 31, 32, 00 (161) 

are linearly independent and span the Lie algebra of U (4) but they give only three disjoint 
sets, viz., {(01), (02), (03)}, {(10), (20), (30)} and {(11), (22), (33)}, containing each 3 
commuting operators, where here again (ab) stands for X a Z h . However, it is not possible 
to partition the set (11611 ) in order to get a decomposition similar to (1 145b . Nevertheless, 
it is possible to find another basis of u(A) which can be partitioned in a way yielding a 
decompostion similar to (1 145b . This can be achieved by working with tensorial products 
of the matrices X a Z b corresponding to p = 2. In this respect, let us consider the product 
u ax bx ® u a2 b 2 , where w aA with i — 1, 2 are Pauli operators for p = 2. Then, by using the 
abbreviation (01610262) for ® u a2 b 2 or X ai Z bl X a2 Z b2 , it can be checked that the 
five disjoint sets 

{(1011), (1101), (0110)} (162) 
{(1110), (1001), (0111)} (163) 
{(1010), (1000), (0010)} (164) 
{(1111), (1100), (0011)} (165) 
{(0101), (0100), (0001)} (166) 
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consist each of 3 commuting unitary operators and that the Lie algebra su(A) is spanned 
by the union of the 5 sets. It is to be emphasized that the 15 operators (|162| )- (|166l) are 
underlaid by the geometry of the generalized quadrangle of order 2 11901 . In this geometry, 
the five sets given by (|162|) - (|166l ) correspond to a spread of this quadrangle, i.e., to a set 
of 5 pairwise skew lines ||90l . 

The considerations of the counterexample can be generalized to d := d\d 2 ■ ■ ■ d e , e 
being an integer greater or equal to 2. Let us define 

uab ■= u aibl © u a2b2 © . . . ® u aebe A := a x , a 2 ,...,a e B :— b±, b 2 ,...,b e (167) 

where w ai 6 1; u a2 b 2 , . . . , u ae b e are generalized Pauli operators corresponding to the dimen- 
sions di, d 2 , • ■ • , d e respectively. In addition, let qi, q 2 , ■ ■ ■ , q e be the g-factor associated 
with d\, d 2 ,...,d e respectively (qj := exp(2iri/dj)). Then, Results 10, 11 and 12 can be 
generalized as follows. 

Result 13. The operators uab are unitary and satisfy the orthogonality relation 



Tr £ ( rflrf2 ... rfe ) 



[uab] U A 'B> 



did 2 . . .d e 5 A ,A> &b,b' (168) 
where 

&A,A' '■= fia 1 ,a' 1 fia 2 ,a' 2 ■ ■ ■ ^a e ,a' e 5 B ,B' '■= ^61,6^62,6^ • • • $b e ,b' e - (169) 

The commutator [uab, ua'b 1 ]- an d the anti-commutator [uab, v>a'B']+ °J u ab and ua'b' 
are given by 

[U A B, U A 'B']t = [Y\. q 3 bjaj ^ II q 3 ^ ) UA " B " ( 17 °) 
\j=l j=l J 

with 

A" := 01 © a[, a 2 © a 2 , . . . , a e © a' e B" := b x © b[, b 2 © b' 2 , . . . , b e © b' e . (171) 

The set {uab '■ A,B E © Za 2 © ... © ^d e } of the d\d\ . . . d\ unitary operators uab 
form a basis for the Lie algebra u(d\d 2 . . . d e ) of the group U (d\d 2 . . . d e ). 

The operators uab may be called generalized Dirac operators since the ordinary Dirac 
operators correspond to specific u ai6l © u a2b2 for d\ — d 2 — 2. 

In the special case where d\ = d 2 = . . . = d e = p with p a prime integer (or 
equivalently d = p e ), we have [uab, ua'b>]~ = if and only if 

e 

"J'j l >j<>'j <• (172) 
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and [uab, ua'B']+ = if and only if 



^a^e^^-p (173) 

3=1 

so that there are vanishing anti-commutators only if p = 2. The commutation relations 
given by (|170|) - (|171| ) can be transcribed in terms of Lagrangian submodules ||65l|9T]|. F° r 
d = p e , there exists a decomposition of the set {uab ■ A,B E Z® e } \ { 1} that corresponds 
to a decomposition of the Lie algebra su(p e ) into p e + 1 abelian subalgebras of dimension 

P e -l[l25l|73l92l|93l|93l95l. 

4.4 Generalized Pauli group 

Let us define the d 3 operators 

Wabc ■= q a x b z c = q a u bc a,b,ce Z d . (174) 
The action of w abc on the Hilbert space £ (d) is described by 

w abc \k) = q a+kc \keb) ^w abc \j,m) = q a+ ^ c \j,m®b). (175) 
The operators w abc are unitary and satisfy 

Tr f (d) [{wabc) ] w a iy c ] = q a '~ a d 5 b ,v 8 c>d (176) 

which gives back (II 141) for a = a' = 0. 

The product of the operators w abc and w a i b > c > reads 

WabcWa'b'd = VJ a "b"c" d" = d @ CL Q cb' b" = b © b' c" — C © c' . (177) 

The set {w abc : a, b, c G Z rf } can be endowed with a group structure. In the detail, we 
have the following 

Result 14. The set {w abc : a, b, c G Z^}, endowed with the internal law < \1 77\) . is a finite 
group of order d 3 . This nonabelian group, noted Yld and called generalized Pauli group in 
d dimensions, is nilpotent (hence solvable) with a nilpotency class equal to 2. The group 
lid is isomorphic to a subgroup ofU(d) for d even or SU (d) for d odd. It has d(d + 1 ) — 1 
conjugacy classes (d classes containing each 1 element and d 2 — 1 classes containing 
each d elements) and d(d + 1) — 1 classes of irreducible representations (d 2 classes of 
dimension 1 and d — 1 classes of dimension d). 
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A faithful three-dimensional representation of 11^ is provided with the application 

A o o\ 

n d -> GL(3, Z d ) : u> afec ^ I b 1 . (178) 




This is reminiscent of the Heisenberg-Weyl group g^|97l 1910193 EXE EDI Indeed, the 
group Hd can be considered as a discretization HW(Zd) of the Heisenberg-Weyl group 
HW(M), a three-parameter Lie group. The Heisenberg-Weyl group HW(M), also called 
the Heisenberg group or Weyl group, is at the root of quantum mechanics. It also plays 
an important role in symplectic geometry. The group U d was discussed by Stovfcek and 
Tolar HI 0211 in connection with quantum mechanics in a discrete space-time, by Balian 
and Itzykson (731 in connection with finite quantum mechanics, by Patera and Zassenhaus 
(741 in connection with gradings of simple Lie algebras of type A n -\, and by Kibler [|25l 
in connection with Weyl pairs and the Heisenberg-Weyl group. Note that the discrete 
version HW(Z) of HWQSL) was used for an analysis of the solutions of the Markoff 
equation [110311 . Recently, the discrete version HW\L p x (Q p /Z p )] was introduced for 
describing (p-adic) quantum systems with positions in 7h v and momenta in Q p /Z p (.104]. 

As far as HW(7id) is concerned, it is to be observed that a Lie algebra n d can be 
associated with the finite group ILj. This can be seen by considering the Frobenius algebra 
of Ud (see HI 051 for the definition of the Lie algebra associated with an arbitray finite 
group). Then, the Lie brackets of %d are 

[W abc , Wa'b'c'}- = Wap-f - W a <f3>y (179) 

with 

a = a®a'ecb' f3 = b®b' 7 = c © d a' = a © cb' be' (3 f = (3 7' = 7. (180) 

The algebra 71^, of dimension d 3 , is not semisimple. It can be decomposed as the direct 
sum 

d 2 d-1 

7r d ~l+|n(l)l+|n(rf) (181) 

1 1 

which contains d 2 Lie algebras isomorphic to u(l) and d — 1 Lie algebras isomorphic to 
u(d). The Lie algebra u(d) spanned by the set {u a b : a, 6 G Z d } is one of the subalgebras 

Of 7T d . 

The group n d (noted P d in (231 ) should not be confused with the Pauli group V n 
on n qubits spanned by n-fold tensor products of ia = <J X and a z used in quan- 
tum information and quantum computation [11061 1107H . The Pauli group V n has 4 n+1 
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elements. It is used as an error group in quantum computing. The normaliser of V n 
in SU(2 n ), known as the Clifford (or Jacobi) group Cli n on n qubits, a group of order 
2n 2 +2n+3 fT™ =1 (4 J — 1), is of great interest in the context of quantum corrector codes 
SH QUI DM Q3D1 ttm [ml El. In addition to Cli n , proper sugroups of Cli n having 
V n as an invariant subgroup are relevant for displaying quantum coherence H114L The 
distinction bewteen V n and 11^ can be clarified by the example below which shows that 
n 2 is not isomorphic to V\. In a parallel way, it can be proved that the groups n 4 and P 2 
(both of order 64) are distinct. 

Example 9: d — 2. The simplest example of n d occurs for d = 2. The group n 2 has 
8 elements (±7, ±x, ±y := ±xz, ±z) and is isomorphic to the dihedral group D A . It 
can be partitioned into 5 conjugation classes ({/}, {—I}, {x, —x}, {y, —y}, {z, —z}) and 
possesses 5 inequivalent irreducible representations (of dimensions 1, 1, 1, 1 and 2). The 
two-dimensional irreducible representation corresponds to 

± 1 \— > ±o" ± x i— > ±a x i y i — > ^fio-y ± z i— > ±a z (182) 

in terms of the Pauli matrices a x with A = 0,x,y,z. The elements e x := x, e 2 ■= y 
and e 3 := 2; of il 2 span the four-dimensional algebra A(l,— 1,0) = Ni, the algebra of 
hyperbolic quaternions (with ef = — e\ — e| = 1 instead of = e 2 = e 3 = — 1 as for 
usual quaternions). This associative and noncommutative algebra is a singular division 
algebra. The algebra Ni turns out to be a particular Cayley-Dickson algebra A(c\, c 2 , c 3 ) 

111 151 . Going back to Il 2 , we see that not all the subgroups of Il 2 are invariant. The group 

11 2 is isomorphic to the group of hyperbolic quaternions rather than to the group Q of 
ordinary quaternions for which all subgroups are invariant (the group Q can be realized 
with the help of the matrices ±a , ±ia x , ±ia y , ±icr z ). Like Q, the group IT 2 is ambivalent 
and simply reducible in the terminology of Wigner fl2]. Indeed, Il 2 is the sole generalized 
Pauli group that is ambivalent. 

To end up with this example, let us examine the connection between n 2 and the 1 
qubit Pauli group V\. The group V\ has 16 elements (±cxa, ±i<J\ with A = 0,x,y, z). 
Obviously, Il 2 is a subgroup of index 2 (necessarily invariant) of V\. The group V\ 
can be considered as a double group of n 2 or Q in the sense that V\ coincides with 
Il 2 [J in 2 = Q |J iQ in terms of sets. Therefore, the group table of V\ easily follows 
from the one of Il 2 or Q. As a result, the numbers of conjugation classes and irreducible 
representation classes are doubled when passing from Il 2 or Q to V\. 
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5 Closing remarks 



Starting from a nonstandard approach to angular momentum and its transcription in terms 
of representation of SU(2), we derived (in an original and unified way) some results about 
unitary operator bases and their connection to unitary groups, Pauli groups and quadratic 
sums. These results (either known or formulated in a new way) shed some light on the 
importance of the polar decomposition of su(2) and cyclic groups for the study of unitary 
operator bases and their relationship with unitary groups and Pauli groups. In particular, 
the latter decomposition and the quadratic discrete Fourier transform introduced in this 
work make it possible to generate in prime dimension a complete set of MUBs given by a 
single formula (Eq. (l20l)). From the point of view of the representation theory, it would be 
interesting to find realisation on the state vectors (l20l) on the sphere S 2 thus establishing 
a contact between the {j 2 , t> 0a } scheme and special functions. 

To close this paper, let us be mention two works dealing with a la Schwinger unitary 
operator bases in an angular momentum scheme. In Ref. HI 161 , unitary operator bases 
and standard (discrete) quantum Fourier transforms in an angular momentum framework 
proved to be useful for spin tunneling. In addition, d-dimensional generalized Pauli ma- 
trices applied to modified Bessel functions were considered in an angular momentum 
approach with j = (d — l)/2 — > oo H117H . 

Appendix A: A polar decomposition of su{2) 

In addition to the operator v ra , a second linear operator is necessary to define a polar 
decomposition of SU (2). Let us introduce the Hermitian operator h through 



3 




ml. 



(183) 



m=-j 



Then, it is a simple matter of calculation to show that the three operators 



3+ ■= hv. 



ra 



3- ■= {v ra ) ] h j z := - [h 2 - {v ra ) ] h 2 v ra \ 



(184) 



satisfy the ladder equations 



j+\j,m) = g +(j - m+s - 1/2) V(j -m)(j+m + m + 1) 



(185) 
(186) 



j_|j,m) = g- (j - m+s+1/2) V(j +m)(j-m + l)\j, m - 1) 



and the eigenvalue equation 



3z\j,m) = m\j,m) 



(187) 
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where s = 1/2. Therefore, the operators j + , j_ and j z satisfy the commutation relations 



+3+ 



\jz,j-)- = -J 



2]z 



(188) 



and thus span the Lie algebra of SU(2). 

The latter result does not depend on the parameters r and a. However, the action 
of j + and j- on \j,m) on the space S(2j + 1) depends on a; the usual Condon and 
Shortley phase convention used in spectroscopy corresponds to a = 0. The writing of 
the ladder operators j + and j_ in terms of h and v ra constitutes a two-parameter polar 
decomposition of the Lie algebra of SU(2). This decomposition is an alternative to the 
polar decompositions obtained independently in If24l II 181 [TT9l 11201 [T2TTI . 

5.1 Appendix B: A quon approach to su(2) 

Following [122], we define a quon algebra or g-deformed oscillator algebra for q a root 
of unity. The three operators a_, a + and iV a such that 



a_a+ - ga+a_ = I [N a , a±]_ = ±a± (a±) k = N% = N a (189) 



define a quon algebra or g-deformed oscillator algebra denoted as A q (a^, a + , N a ). The 
operators a_ and a + are referred to as quon operators. The operators a_, a+ and A^ a are 
called annihilation, creation and number operators, respectively. 

Let us consider two commuting quon algebras A q (a-, a + , N a ) = A q (a) with a = x,y 
corresponding to the same value of the deformation parameter q. Their generators satisfy 
equations (11891) and (11901 ) with a = x,y and [X, Y]_ = for any X in A q (x) and 
any Y in A q (y). Then, let us look for Hilbertian representations of A q (x) and A q (y) 
on the A;-dimensional Hilbert spaces F{x) and F(y) spanned by the orthonormal bases 
{\n\) : Hi = 0, 1, . . . , k — 1} and {\n 2 ) : n 2 = 0, 1, . . . , k — 1}, respectively. We easily 
obtain the representations defined by 



where 





(190) 



x + \ni) = \ni + 1) : 
a;_|m) = [m] g |ni - 1) 
iV x |ni) = ni|ni) 



x+ \k - 1) = 
) x_|0) = 



(191) 
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and 



y+\n 2 ) = [n 2 + 1] q \n 2 + 1) - 1) = 

y_K) = K-l) J/-|0) = 

iV> 2 )=n 2 |n 2 ) (192) 

for A q (x) and respectively. 

The cornerstone of this approach is to define the two linear operators 

h := ^N x (N y + l) v ra := s x s y (193) 

with 

Sx : = ^(W^+e^ 1 . (x,)^ 1 (194) 

Sy := y _g-^--^)/2 + e *«r/2^^( y+ )*-l (195) 

F ~~ X J c/ ! 

where 

oGR r = tt(A; - l)r r G R. (196) 
The operators h and v ra act on the states 

\ni,n 2 ) :— \ni) ® |n 2 ) (197) 

of the /c 2 -dimensional space jF fc := ^"(x) ® ^(y). It is immediate to show that the action 
of h and v ra on jF fc is given by 

h\n 1: n 2 ) = Vwi(w2 + 1)K,"2) ^ = 0, 1, 2, • • • , fc - 1 i = l,2 (198) 

and 

v ra \n 1 ,n 2 )=q an *\n 1 + l,n 2 -l) n^k-l n 2 ^ (199) 

^a|A;-l,n 2 ) = e^ / V a( ^ 1 ~ ri2)/2 |0,n 2 -l) n 2 ^ (200) 

«ra|ni,0) = e^ / y (fc+ni)/2 |ni + l,A;-l) m ^ fc - 1 (201) 

Vra\k- 1,0) =e**-|0,ib- 1). (202) 
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The operators h and v ra satisfy interesting properties: the operator h is Hermitian and the 
operator v ra is unitary. 

We now adapt the trick used by Schwinger [ 12311 in his approach to angular momen- 
tum via a coupled pair of harmonic oscillators. This can be done by introducing two new 
quantum numbers J and M defined by 

J := -(m +n 2 ) M := i(m -n 2 ) \ JM) := \ J + M, J - M) = \m, n 2 ). (203) 
Note that 

j:=\{k-l) (204) 

is an admissible value for J. Then, let us consider the /c-dimensional subspace e(j) of the 
& 2 -dimensional space J-{x) ® spanned by the basis { \j, m) : m = j, j — 1, . . . , — j}. 
We guess that e(j) is a space of constant angular momentum j. As a matter of fact, we 
can check that e(j) is stable under h and v ra . In fact, the action of the operators h and v ra 
on the subspace e(j) of jF fc can be described by 

h\j, m) = yj {] + m)(j-m + 1) \j, m) (205) 

and 

v ra \j, m) = S m ,e^ r \j, -j) + (1 - 5 mJ ) q^ m)a \j, m+1) (206) 
in agreement with Eq. (11831) and with the master equation ©. 
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